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ABSTRACT 


In salvo firing, the smallest ballistic dispersion is 
not always most desirable. Deliberate increase of the 
ballistic dispersion can improve the probability of des- 
troying the target. Our concern in this dissertation is 
mie Obtimization of such “artificial” dispersion in two- 
dimensional salvo models. In some cases no closed form 
solution is available, but we are able to offer efficient 
methods for the computation or approximation of the salvo 
kill probability. In other cases we are able to derive 


approximate formulae for the optimal ballistic dispersion. 
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Let us begin with a classical example presented in Chap- 
ter 6 of the textbook by Morse and Kimball [Ref. 1]. Suppose 
an airplane carries two bombs to attack a railroad track. 

The plane flies along a course perpendicular to the track, 
and drops the bombs. There are three ways to drop two 

bombs: (1) together in salvo, (2) spaced a suitable distance 
apart aiming the midpoint of the pattern at the center of 

the track, and (3) dropping each on a separate run over the 
eacget. 

We intend to illustrate the method of calculating the 
probability of track destruction in the three cases, the 
purpose being to determine which method of attack is the 
best. 

Let us consider Case (1), salvo bombing. The two bombs 
leaving the plane simultaneously will hit the ground at some 
distance apart; the impact point being random. For the time 
being the random vector of the impact position from the 
center of impact is termed the ballistic error. The aiming 
is also not free from error. The plane aims at a point on 
the center of the track, but the center of impact will 
deviate from the aim point. This deviation, or aiming error, 
is also random. The impact point of each bomb is, there- 


fore, composed of two errors; first the aiming error which 





1s common to both bombs, and second the ballistic error 
which varies from one bomb to another. 

If the variance of the ballistic error is small while 
that of the aiming error is large, the two bombs will land 
almost at the same point. Therefore it is expected that 
either both will hit or neither will hit the track. To avoid 
the latter situation and to improve the probability of des- 
truction, it is better to spread the landing pattern of the 
bombs. One way is the method of pattern bombing described 
as Case (2), another is salvo bombing employing bombs 
with a suitable value of variance of the ballistic error. 
Our Main concern in this dissertation is maximization of the 
probability of salvo destruction by choosing a suitable 
value of the variance of the ballistic error. In the follow- 
ing we call it the optimal ballistic dispersion. 

The problem described above is not a new one. In an 
interesting talk before a conference at the Ballistic Re- 
search Laboratories held in March 1955 [Ref. 2], Merritt 
and King stated that an approximate formula for the optimal 
ballistic dispersion was derived by two Englishmen as early 
as in 1936. Following the work in WW II, a number of arti- 
cles were published on Boch the salvo and pattern firing 
models, but since the early 1970's, it seems that these 
models have attracted less attention. Further detail can be 
found in two excellent review papers by Eckler [Ref. 3] and 


Eckler and Burr [Ref. 4]. However, little is known as to the 





optimization of the ballistic dispersion in salvo models. 
It is the calculation or approximation of the optimal 


ballistic dispersion that is the subject of this dissertation. 


A. SALVO FIRING 
In the following chapters we deal with only two-dimensional 
salvo models. The reason we adopt the two-dimensional model 
is that it is the most frequent case and plays an important 
role in real world applications. One- or three-dimensional 
theory can also be developed to parallel our investigation. 
Suppose that there is a target in a two-dimensional space, 
and a salvo of n weapons is delivered against the target. 
Delivery error relative to the target is assumed to be com- 
posed of two parts, the aiming error and the ballistic error. 
First we aim at the target. With this aiming, the center 
of impact point of n weapons is determined. Let us adopt a 
Cartesian coordinate system (x,y) such that the origin coin- 
cides with the center of impact. The aiming cannot be per- 
fect. Let the components of the aiming error be -U and -V. 
Then the position of the target with respect to our coordinate 
system 1s given by (U,V). It is assumed that the joint 
probability density of U and V exists and is given by £, (u,v). 
Now, a Salvo of n weapons is fired against the target 
after (U,V) takes a value (u,v) which is unknown to us and 
only is predictable in a probabilistic sense. The components 
of the impact point of the ith weapon are denoted as Xs and 


Y.- It is assumed that (X.Y) are independent and identically 





distributed random variables with joint density £,(x,y) . 
mtg. |, a typical geometry is shown. The target is at 
(u,v), and the impact points of four weapons are scattered 
mareund the origin. | 

As to the target we adopt the so-called point target 
concept: It 1s assumed that the target is either completely 
destroyed or else undamaged by each weapon. We neglect any 
possible partial damage and its cumulative effect. The 
probability that a weapon landing at (x,y) destroys the tar- 
get at (u,v) is a function of u-x and v-y, denoted as 
D(u-x,v-y) and is called the damage function. 

The probability of destroying the target is calculated 
in two steps. First, suppose that the random variables 
U and V take some values u and v. The conditional proba- 
bility that the ith weapon destroys the target given U = v 


amc VY = Vv is 


p(u,v) = J f D(u-x,v-y) £, (x,y) dxdy : Ciel): 


The impact points Xx. and v5 are all assumed to be independent, 
so the conditional survival probability of the target given 
Me= uand V= vis {1 - p(u,v)}>. Therefore, we obtain the 
probability that the target is destroyed by the salvo of n 
weapons--salvo kill probability--by averaging the conditional 


survival probability with respect to the distribution of u 








lanl eibuat=s Ae 


10 





ev, and then subtracting it from unity: 


oO oO 


me i-f f ti- p(u,v)} £, (u,v) dudv Us 


=e OO —_ 


B. UPPER BOUND ON THE SALVO KILL PROBABILITY 

In the example mentioned earlier, Morse and Kimball ob- 
served that the probability of destroying the track by drop- 
ping two bombs in salvo is always smaller than the destruc- 
tion probability of Case (3), namely the case in which 
the two bombs are dropped independently on separate passes. 
Similar observations are also pointed out in the review paper 
by Eckler and Burr. The observation is valuable, but there 
has been no proof given to this fundamental property dis- 
covered numerically. Indeed, we have the following proposi- 
tion without any specific assumptions on f, (u,v), £. (x,y) 


ema D(x,y) . 


Peoposition 1.1 
Let the salvo kill probability with n weapons be denoted 


as Ph 


OO 


= 1 - fe f a = p(u,v) }'£, (u,v) dudv aes) 


Then, 


ded. 





Peele (= ps)” 


a4 1 (igs) 


Namely, salvo Ee 1S always smaller than the kill probability 


of n independent tries. 


Proof 

The proof is immediate from a theorem in the book by 
Hardy and others [Ref. 5]. Define a mean M. of a non-negative 
function g(x) for a positive r as, 


M.(g) = (fog dx}l/¥ , 


where £ = £(x) is a weighting function, positive everywhere 


and, 
fied = 


Theorem 192 from Ref. 5 tells us that if 0 < r < § and 


M,(g) is finite, then 


unless g(x) 1S a constant. 
The weighting function £ corresponds to f, (u,v) i Out 
paeposition, and g corresponds to our 1-p(u,v). Therefore, 


for n > 1, we have 


We 





M, (1-p) < M_ (1-p) 
or 

(M,(1-p)}" < {mM (1-p)}" , 
which means 


ae ee ee 


(1-P, we Ui 


Here is an outline of the succeeding chapters: In Chap- 
ter II, we investigate the simplest salvo model which is 
characterized by circular normal errors and circular Gaussian 
damage function. After that, an approximate formula for the 
optimal ballistic dispersion is introduced in Chapter III; 
we regard the material in Chapter III as the "core" of this 
dissertation. In Chapter IV, the approximate formula is 
applied to the salvo model of Chapter II, and its accuracy 
is studied. In Chapter V, an approximate method of calcu- 
lating the salvo kill probability of the salvo model with 
the general damage function is presented; this case is more 
difficult than the one where the damage function is circular 
normal, so we resort to further approximations. In Chapter 
VI, we investigate the salvo model with the so-called cookie 
cutter damage function, relying heavily on the results of 
Chapters III and V. Finally, salvo model without circular 


symmetry is discussed in Chapter VII. Our concern throughout 


is 








memeene Optimization of artificial dispersion. In some cases 
no closed form solution is available, so we offer efficient 
fecnoas for the computation or approximation of salvo kill 

probability with a given dispersion. In other cases we are 


able to offer closed form solutions. 
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heweolveies SALVO MODEL 


In this chapter we deal with the simplest salvo model 
with circular normal errors and a circularly symmetric 
Gaussian damage function. 

Let us adopt a Cartesian coordinate system (x,y) the 
Origin of which coincides with the center of impact of the 
mn weapons. The target position with respect to this coor- 
dinate system is denoted as U and V in the x and y direction, 
respectively. (U,V) is a random bias common to all the n 
weapons. In this chapter (U,V) is assumed as a circular normal 
variate; namely, we assume that U and V are independent and 
identically distributed normal random variables with mean 0 


and variance a. The joint density of U and V is, therefore 


e d @)) 





f, (u,v) = 


The impact point of the ith weapon is denoted as (X.,Y.). 
It is assumed that (X.,¥5), 1 =1,2,...,n are independent and 
identically distributed circular normal random variables 


with joint density function 





£., (x,y) = ze : (22) 


dee 





As to the damage function, we assume that the so-called 
Gaussian damage function with circular symmetry: If the 
target is at (u,v) and a weapon impacts at (x,y), then the 


destruction probability is given by 


_ (uex) 74 (vey)? 
Der 
D(u-x,v-y) = e : (223) 


In the following, let us call this model the Simple Salvo 
Model. 

According to Grubbs [Ref. 6], the Simple Salvo Model was 
first studied as early as in 1953 by H. K. Weiss in BRL Re- 
port No. 879, "Methods for Computing the Effectiveness of 
Area Weapons." The formula for the salvo kill probability 
(2.7) given in the next section is attributed to him. Later, 
Breaux [Ref. 7] found that this formula was not suitable for 
computation when n is large, and gave another method of com- 
putation. However, it seems that there is no published re- 
search investigating this model in full depth; the studies 
usually come to an end when expressions for the kill proba- 
bility are derived. In this chapter, we investigate first 
the property of the kill probability as a function of param- 
eters involved in the model. Methods for computing the kill 
probability and bounds for the kill probability are also 


dealt with in the later sections. 
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eee VvO KLE PROBABILITY AS A FUNCTION OF n, oO fo} AND a 


1 ad 
The salvo kill probability is derived in two steps as 

is stated in Chapter I. First the conditional kill proba- 

bility by the ith weapon is calculated given that the random 


bias (U,V) takes a value (u,v). 
p(u,v) = ff D(u-x, v-y) £, (x,y) dxdy . Gig) 


The salvo kill probability is then obtained by averaging 


[ert -p(u,v)}" over the distribution of (U,V), as 


co 


eee 2-f f {i-plu,v)}"f, (u,v) dudv . (1.2) 


=O 


Using the assumptions (2.2) and (2.3), we easily calcu- 


late the conditional kill probability (1.1) 





_ (u-x) 7+#(v-y)*% _ x*+y? 
ho a Ber 204 
aay) | = ‘~oe.) 2 dxdy 
2105 —co —co 
and get 
_ a 
2 maar =) 
a Z 
Gulia) = a) e : (2.4) 
Q- FO 


129; 





Bmeseleucing (2-1) and (2.4) into (1.2), the salvo kill 


Peobability is 











_ ur+v" = u-ty 
1 co G0 Ae 2 (a*+0%) n 204 
P = 1- 5>j | t1-—3=—7e toe dudv 
104 —< 800 OQ +o, 


Conversion from Cartesian coordinates into polar coordinates 


gives 
2 2 
2 r Bis je 
l oo . 2 +05) n 20% 
ae) —— C } e r= lig 
0) 0 a +0 
i 2 
From this we get 
Bs a Ae (2.5) 
a P 
where 
a Ze? _ z 2 2 
KH = @ [Os Sas (Ce +05) /01 ; CORRS) 


There are 4 parameters in our model, but the salvo kill 
probability is determined by three, n, O,/a, and C5/o, or 
n, A, and p aS is seen in (2.5). 

If the integrand in (2.5) is expanded in a binomial 


series and integrated term by term, we obtain the Weiss 


le: 





femmuLa for the salvo kill probability. 
CRor =e en) 


P is obviously an increasing function of n, and is de- 
creasing in oat ea The latter is easily verified by in- 


2 


specting (2.5), where \/op is constant for fixed oS /a The 


o$/a”-dependence HOemNOWe Ver, Not Obvious. in (2.5), A is 
constant for a fixed value of of/a*, and there are two factors 
which involve p in the integrand, with p tending to increase 
P in one and to decrease P in the other. 

Table 1 and Table 2 are given here to show the of /a*- and 
04 /a”-dependence of the salvo kill probability, where n is 
kept constant, n = 2 for Table 1 and n = 16 for Table 2. For 
fixed n and oS /ar, P is decreasing in o¢/a as mentioned 
above. However, the o$/a°-dependence is a little different. 
It is observed that, generally, P is decreasing in a3 /a* 1 eng 
small values of a2 /a", but for a larger af /a", P increases 
when aS /a" is increased, reaches a maximum value, and then 
decreases. It is also noted that from the viewpoint of real 
world applications, the case of large o¢/a* with large n is 
especially interesting; e.g., as seen in Table 2, the kill 
probability for n = 16 and o*/a* = 8.0 is improved to 0.5881 
when oS /a* is increased to a suitable value. The maximum P 


is more than 1.7 times the original value for a3 /a = Q. 
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Table 2 Salvo Kill Probability 


n= 16 
De 2 
0,/a 
a3 /a* 5.20 5.5 1.80 2.8 4.09 8.20 
G.ad 1.9990 B.9935 G.9412 G.7335 3.5511 8.3370 
3.25 1.8082 B.9956 G.957G 6.8234 0.6013 0.3784 
B.5d 1.9990 B.9965 B.9659 9.8504 3.6399 6.4128 
G.75 1.8080 3.9968 G.971G 0.8692 B.6702 G.4417 
1.92 1.3990 B.9967 g.9739 G.8824 8.6938 G.4662 
ie 1.25 g.9999 B.9962 3.9752 9.8915 G.7126 9.4871 
| 1.52 G.9997 g.9954 B.9753 3.3978 9.7276 B.5849 
ie. 75 3.9993 G.9941 G.9745 G.9817 0.7394 9.5201 
2.86 B.9985 B.9922 G.9728 B.9838 G.7485 B.533¢0 
| 2.25 B.9972 3.9899 8.9704 8.9845 0.7555 3.5442 
2.5@ G.9954 B.9869 G.9672 8.9838 9.7606 8.5534 
Pe S B.9930 0.9833 3.9635 9.9821 0.7642 0.5612 
3.90 G.990¢G G.9792 g.9591 6.8996 9.7665 8.5678 
3.25 B.9863 @.9745 G.954] B.8962 G.7676 G.5732 
3.52 9.9821 8.9692 G.9486 9.8921 Q.7677 G.5776 
8.75 Q.9772 0.9635 Q.9427 3.8875 B.767@ 9.58190 
«4.868 G.9719 G.9573 B.9364 B.8824 S.7655 G.5837 
4.25 G.9662 G.9587 G.9298 3.8769 0.7634 3.5857 
4.58 0.9597 B.9438 G.9228 G.8712 G.7607 G.587@G 
4.75 G.953¢d 8.9366 B.9155 9.8648 B.7576 3.5878 
5.00 B.9459 G.9291 g.9081 6.8583 G.7541 8.5881 
5.25 B.9386 G.9214 B.9804 6.8516 6.7502 3.588¢ 
5.50 B.9309 9.9135 B.8926 9.8448 6.7462 @.5874 
5.75 9.9231 9.9055 0.8847 9.8378 9.7416 3.5865 
6.38 G.9151 9.8973 8.8767 0.33087 G.7369 8.5853 
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Experience shows that there is one and only one local maximum 
of P as a function of oS /a* fom 2ixea nh and af/a”, but we 
have not succeeded in obtaining a proof. 

The salvo kill probability can be increased by increasing 


2 


03/0 from zero for given n and oS /a if the condition below 


is satisfied. 


Eeoposition 2.1 


The salvo kill probability P for fixed n, 0, anda lis 


1 
not maximal at 2 CMirteancd Only sf 
iL il a 

Tae oe i ers > Jy (2.8) 

where 
2 
XH = a /O4 
Proof 
We will show that 
dP/do | > 0 (25:93) 
p=A 

memeaquivalent to (2.8). Using (2.7), (2.9) is written as 


“t/dyn-L-t/r-tae oon f (ree tA Phe te) at 


n ii (l-e 
0 0 


(2.10) 


ta J 





Now define 


ie 


meee) nf (hoe 74) te f/A-Xt 
0 


Maem the inequality (2.10) is rewritten as 
1121 1) (Ie <a ele Gl) Te a aaa CZlek) 


where I'(1) means the value of the derivative dI/dx at x = l. 


Mee s = exp(-t/A) in the integral, then I(x) is 





AX n-L 
Mmejee= nA {| s ~(1=s) ds 
0 
a SI: sg ces 
7 AX+L Ax+2 A\X+n 
Therefore, 
Cail 3s) re ; am 
Te) 4=1 IP as 


The inequality (2.8) is readily obtained from this 


Saidation with x = 1 and (2.11). Oj 


Proposition 2.2 
A sufficient condition that the salvo kill probability 


eetor fixed n, Cis and a@ is not maximal at 5 =O 1.S 
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i) 
i) 
IA 


o! /O+ n/(e-l) - 1. (Coe) 
Preoot 


pamce 1/(x+tA) 1S decreasing in x, 


qsPul n+l 
Co f 


ine: 1 


n+1+i 
1+ 





io dle 
—— > — dx = In 
+A + 
1 J 5 ss 


Il ~35 
=> 
”~ 
+. 
> 
l 


it5 


ieecthe Fright-hand-side is equal to or greater than unity, 
i.e., lt+n/(1+A) 2 e, then the condition (2.8) in Proposition 
2.1 is satisfied. The relation (2.12) is immediate from the 


above-mentioned inequality. a 


Proposition 2.3 
Assume that (2.8) holds. A necessary condition that the 


salvo kill probability be maximum at some point o., > 0 for 


Ps 

fixed n, 07 and qa 1s 

a Jol joy op eo eee 

P(e)" Gy g(r One (ee i) 

j=1 J P (aay)! 
where 

Se. _ Ce ae: 

Proof 


Differentiate (2.7) with respect to op, and equat it to 


O. The equation (2.13) is immediate. | C 
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Bxample. Case n = 2 
Equation (2.13) is generally too complicated to be solved 
Or 05/0 in closed form. Here we solve it for the simplest 
case, n = 2. 
The inequality (2.10) with n = 2 is 
1 
i 


1p 
1+} ite ue 1, 


and after a bit of algebra, this reduces to 


We cay a= 27 VS -1 4) (2.14) 
Now, assume that ot /a" is large and satisfies (2.14). From 
13) 
2 —2. - 24 oF) = 0 
(o+1) P (p+2) 
Therefore, 
Doe ( Aye 
O5/a = of/i-l = ——, - = 
oe (ese } 
(5) 
ae Goss = 
01/0 = aN 5) 
Sle.) 


This is a parametric expression of maximizing oS /a" as a 
Banmetion of os /a° the curve of which is shown in Fig. 2. The 


value of the parameter p which corresponds to the point 
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on / a e 
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eye = (v5 +1)/2, a5 /a* = 0) is p = 09, where 
op = (v5 = 1)/2 


When op approaches Zero, 


2 


oF /a* ~ Wiles ; oS /o e/a, 


and therefore, the asymptotic form of the curve is 


ee 
J5/a = 5 o,/a : (6221.6) 


B. COMPUTING THE SALVO KILL PROBABILITY 

In the preceding section, the Weiss formula for the salvo 
kill probability was given in (2.7). According to Breaux 
(Ref. 7], this formula is not suitable for calculation for a 
large n because it 1s an alternating series. Breaux pointed 
out that the salvo kill probability could be expressed in 
terms of an Incomplete Beta Function, and that for large n, 
Tang's method [Ref. 8] would be advisable for calculating an 
Incomplete Beta Function. 

First, we present a recurrence formula for calculating the 
salvo kill probability which is, in principle, the same as 
that of Tang's method. Its generalization will be given in 


a later chapter. 
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Proposition 2.4 
Let us denote the salvo miss probability by n weapons as 


O(n) = 1={P. 


ieee (mate 7h) 6° dt . (eee | 
0 


Then Q(n) satisfies the following recurrence equations: 


= n 0 n = 
Q (n) = nto Q(n-1l) + Ato (tT) . Time De cus < CZs) 
OCO 7 = 1, (2ie 19) 

where 

me = rA/po = a°/(a* +04) - 0 = (a* +04) /o4 
Proof 

Erom (2.1/7) 
Q(n) = f meee)" ea at 


foe) 


f Gene 72)" te “at = ay Gene ©! Oy ee ea 
0 0 
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The first term in the right hand side is Q(n-l). In the 


second term, we integrate it by parts, and have 


Qin) = Q(n-1) - 2(1-pe*/PyMeM®]  - Bf (a -pe™*/P) Me™ at: 
0 0 


Q(n-1) + 2(1-y)" - 2 Qin) 


The equation (2.18) is immediate. 0 


It is worthwhile to note that the method presented in 
Proposition 2.4 gives not only the value of Q(n) for a 


specified n but also gives all the Q(n)'s up to Ny. It 


0! 
1S an important advantage of this method over the Weiss 
formula because we are often interested in the salvo kill 
probabilities for several n's. As is shown in Table 3, the 
computational time for Q(n) by our recurrence algorithm is 
about twice the time required for computation of a single 
Q(n) using the Weiss formula. When the salvo kill probabili- 
ties for more than a single n are needed, therefore, use of 


the recurrence algorithm given in Proposition 2.4 is 


recommended. 


We have another efficient algorithm for computing the 


Paevo Kill probability, which 1s good even for large n. 
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Table 3. 





Comparison of the Processing Time for 


Computing the Salvo Kill Probability 


Mean time (microsecond), 


pe = ana et — 055. 


n Weiss formula Recurrence formula 
LO 92 neva 
20 189 383 
30 300 sie) 


Mean time (microsecond), 


=p anden — 20. 


u Weiss formula hecubnencoecOonmura 
O21 90 196 
92 175 
0.9 100 186 


Computation is on the IBM 3033 installed in the Naval 
Postgraduate School. Each value is the mean of 4 trials. 
One trial is 1000 repetitions of calculation, and the 


elapsed time units are divided by 1000 and multiplied by 


26 (microseconds). 


Sv 





Ppeposition 2.5 


The salvo miss probability Q(n) by n weapons 
One ee 9°)" e - at e17) 
0 
satisfies the following relationship 
a n n-k k 
Co) a) ue q(k) , (2.20) 
k=0 


where gq(k) is the salvo miss probability by k weapons with 


p=, 
me ei G@-e0 °) en” at 
0 
Proof 
ey = Of (2 -u) + w(-e 97") 3" 2 at 
0 


al ore) 
) ©) (l-u)” ; u* - (1 -pe Eau ea at 
=0 


Si 








Souollary 2.6 
imessa!vo miss probability (2.17) can be calculated by 


the formula 


nN 
Ot) A, , 19) os 
k=0 
where 
eee = (n= Kk+1)/(1-n) (k+p) , k = 1,2,... 
(2a 2 2)) 
Ay = (l=) 
Proof 


As can be seen from (2.18) and (2.19) with yp = 1, q(Q0) 


is unity, and q(k)/q(k-1) is k/(k+p). Utilizing the fact 
that (7 (Gog = (n-k+1)/k, the corollary can be obtained by 
Supstitution. 0 


Note that all the terms in the expression (2.24) are 
positive, so there is no problem arising from cancellation 
of terms with alternating signs as was observed in the Weiss 


Bormula (2.7). 


In checking the computation, the following P's for special 
cases would be of use. 
(l) The case a, = OQ corresponds to 0 + », and we have 


2 ae eee eet 
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(2) The case O5 = 0 corresponds to nu = 1, 9 =A, and so 


from the recurrence formula, we have 


ighivige=It)) OOo If 


SOG srsimese lal) (2.24) 


(3) The case a“ +05 = cae 


is easily integrated to give 


Tae cLormuta “(225)' with o = 1 


= a ae 7 > n+l 
pS Te See Wl) Sie (2.25) 
where 
uo= a*/(a* 4 oS) , A = a“ /os 


See SOUNDS TO THE SALVO KILL PROBABILITY 
In Chapter I it was pointed out that an upper bound of 
the salvo kill probability is given by the kill probability 


of n repeated independent shots. 


where 


oO 


P = 1-([f f {1 - p(u,v) }£, (u,v) dudv]” 


we OO 


Under the assumptions of the present chapter, we have 
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SE 


Proposition 2.7 


An upper bound to the salvo kill probability is given by 


Ree 
= (l= ee (2.26) 


Ul 
Tl 


where 
= Zz Zs 2 2 2 


As to a lower bound, we have 


Peeposition 2.8 (Merritt and King) 


A lower bound to the salvo kill probability is given by 
= - Pp 
P = 1 = p (ny) Bo r (2e27)) 
where te) is the Incomplete Gamma Function, 


Bre 
f = ,e7t ae. , 


P(e) = 

and 
2 <a 7 7 
wo = a /(a° +05) , i OO, 


Table 4 and Table 5 show how the true kill probabilities 


are bracketed by these two bounds. The triplets of entries 
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Table 4 


9.4898 
2.4898 
9.4353 


9.4375 
9.4375 
@. 3935 


e951 
Broo 5 1 
‘DS SUS sts! 


9.3600 
0.3600 
e209) 


Salvo Kill Probability and Its Bounds 


0.50 


0.7506 
eT See 
0.6187 


0.6409 
8.6333 
0.5443 


9.5556 
6.5524 


9.4898 
9.4881 


9.4375 
J.4365 
G25922 


6.3951 
9.3944 
0.3589 


9.36000 
G23 596 
Oe3291 


0.3306 
G3 503 
0.3045 


rn = 


Dies, 
0,/o 


1.00 


3.6095 
6.5216 


“a Pie Bo Bolo) 
6.5417 
9.4715 


6.4825 
6.4267 


0.4375 
9.4333 
0.3882 


O23951 
O.3925 


0.3600 
Jiesoo3 
6.3274 


9.3306 
9.3294 
0.3032 


0.3056 
6.3048 
0.2823 


Z 


2.00 


9.4898 
0.4562 
Ooo Ey 


9.4375 
9.4167 
J200/9 


Bas Io 
BSS S\P4S) 
6.3429 


9.3600 
0.3524 
‘D) Auge) oP 


6.3306 
SSW S) 
reo) 7 


O.eb2) 
6.2784 


9.2840 
Oe 2eL5 
0.2611 


Diao 3 
9.2635 
9.2457 


ps 
Q 
© 


ho 
186) 
J 
UJ) 


Q AAQ WAA RAA AAA @f2QA AAA AQABA 
th 
GO 
78) 
00 


8.00 


9.1994 
9.1724 
O2lsi27 


3.1900 
Orebe 22 
et 522 


G.1814 
@.1689 
9.1542 


D136 
9.1643 
Oe 56 


0.1664 
Orne 
9.1482 


Ceo 
6.1542 
9.1444 


Seas) Eis. 
6.1492 
0.1464 


0.1479 
9.1443 
9.1365 


Triplets are, from the top, P, the upper bound, 


P, 
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the true value, and P, 


the lower bound. 





Table 5 


0.00 


1.8000 
1.0000 
997 


me 2985 
0.9885 
22] 2 


G.9900 
6.9900 
2.9517 


9719 
eeo7 19 
me 2592 


6.9459 
9.9459 
229305 


e915] 
ell 
O25 9683 


6.8819 
@.8819 
9.8647 


6.8481 
6.8481 
8.8310 


Triplets are, 


P, the true value, 


Salvo Kill Probability and Its Bounds 


0.508 


939997 
0.9967 
0.9944 


0.9954 
0.9922 
0.9862 


Oe 9321 
O.9792 
0.9685 


he SMB) 3) 7 
B.29573 
0.9432 


OB. 9309 
Ooo 
0.9130 


0.8987 
0258973 
6.8802 


Q@.8656 
078039 
0.8468 


0.8313 
0.8305 
9.8136 


n= 


a8 /a 


1.80 


O22I9G0 
O27 22 
0.9688 


0.9908 
6.9728 
9.9644 


So Seite, 
Bio Do 1 
0.9469 


0.9459 
0.93604 
Oioe > 


Oo 151 
0.9081 
6.8917 


0.8819 
0.8767 
0.8597 


6.8481 
9.8441 
Oe 27 3 


0.8147 
Betsy baby 
O27 953 


IES) 
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P, 


2-00 


9.99900 
0.8824 
0.8750 


9.9719 a 


9.9038 
058935 


@.9459 
O02 5996 
0.8804 


Ge shinee 
0.8824 
9.8074 


Oe Gel 
5) 3 (eras) 6 
0.8424 


9.8481 
6.8307 
@.8145 


@.8147 
0.8016 
9.7857 


9.7824 
@.7724 
9.7576 


from the top, je, 
and 


0.0955 
9.6867 


Oe ow 
@.7485 
0.7388 


0.8819 
9.7665 
9.7546 


0.8481 
0.7655 
Dot oZ22 


@.8147 
0.7541 
0.7462 


0.7824 
9.7369 
Oei2 29 


OB. 7515 
Oirey Lod 
0.7030 


Oa i222 
9.6952 
0.6819 


8.00 


6.8147 
0.46062 
@.4611 


0.7824 
Geos 
a8 Saye: 


4.6 els 
0.5678 
G2 5590 


Ons 222 
Ge Byes) Uf 
O50 39 


9.0945 
5.6 Sele lh 
Ga) UIs 


0.6684 
Be5G55 
0.5749 


9.6439 
9.5786 
O62. 5677 


0.6209 
0.5680 
See) Te, 


the upper bound, 
the lower bound. 





in the table are, from the above, the upper bound P, true 

feeeeo P, and P, the lower bound. It is observed that 

P-P is not small for smaller values of sa So) but gets smaller 
when 05 /a* is increased. If a4 /a* is small, P and P are 

paigce. | On the other hand, P is fairly close to P over wide 
ranges of o¢/a° and anja Figures 3 and 4 illustrate this. 
memcectis reasonable that Merritt and King used P as an 
approximate formula for P in their study on the optimal 
ballistic dispersion, because computation was really a 

problem at that time. Nowadays, however, direct computation 


Sar is much more natural. 


By) 





1.00 


0.98 


One Its: 


0.94 


1.0 


2.0 


Figure 3, 
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Figure U4. 
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III. OPTIMIZATION OF SALVO FIRING 


In Chapter II we have learned the following: The salvo 


kill probability is a decreasing function of oy if other 


parameters are kept constant. Therefore, to obtain higher 


salvo kill probability, it is necessary to make J, as small 


as possible. The salvo kill probability as a function of 


C5 is, however, complicated. If O71 is sufficiently small, 

the salvo kill probability is monotone decreasing in Tos 

Buc if o1 is not so small and satisfies the condition (2.8), 
then the salvo kill probability is not maximal at S54 = Oe 

but takes a maximum value at some J4 Owen “OrMer words, tHe 


smallest ballistic dispersion is not always most desirable. 
Deliberate increase of the ballistic dispersion is encoun- 
tered in shotguns and in related military weapons. In this 
chapter we present an approximate formula for the optimal 
ballistic dispersion. We will first briefly sketch Walsh's 
theory [Ref. 9] since our approximate formula is based on 
it. In Section B our approximate formula is presented. 


Its accuracy will be studied in Chapter IV. 


A. WALSH THEORY 

The Walsh model is based on assumptions similar to those 
of Chapter I. A salvo of n weapons is fired at a target. 
The impact point error is composed of two parts, a random 


bias which is common to all n weapons, and a round-to~round 
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error. A Cartesian coordinate system is chosen so that its 
Origin coincides with the center of impact of n weapons. 
The coordinates of the random position of the target are 
denoted as U and V with respect to this system. The joint 
probability density of U and V is given by f, (u,v). The 
impact points of weapons (X5,Y;) are independent and iden- 

| tically distributed random variables with joint density 
£, (x,y). It is assumed that impact of a weapon at (x,y), 

| given the target is at (u,v), destroys the target with 

| probability D(u-~-x,v-y) independently of other weapons. 


The lethal area of a weapon is denoted as A: 
J f D(x,y)dxdy = A 


The salvo kill probability is given by (1.2), 


co 


P = l1- f f a - p(u,v) }' £, (u,v) dudv , G2) 


where p(u,v) is the conditional destruction probability by 


a single weapon given that the target is at (u,v). 


oO co 


p(u,v) = f  f D(u-x,v-y) f(x,y) dxdy . (si) 


Note that p(u,v) can be integrated to give A. 
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ff pluyvidudv = af f £,(x,y)axdy = A. (3.1) 


The problem is to derive the aye Rapa ballistic disper- 
sion, ae Instead of solving the original problem, Walsh 
solved a revised problem that was different from the original 
one in that the variable is not af itself but rather the 


meron p(u,v) . 
Proposition 3.1 (Walsh) 


foe) foe) 


max eS = | = f i) {1 - p(u,v) }'£, (u,v) dudy (3.2) 
p(--) —CcO =Cco 
subject to 
le erpltu,v)dudy = A; 
(33) 
Sina) Sau 
has a solution 
i 
= é n-l 
ol (etuv)’ (u,v) e E 
p*(u,v) = (3.4). 
0 (u,v) « E 
where 
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m= { (u,v) |£, (u,v) os aan (35) 
C is positive and is given by 
f f pay) dudy 8=" A. (3.6) 
E 
The corresponding salvo kill probability is 
eer i) to, (uv) sel dudy - (3.7) 
E 


In case the random target bias iS given by a normal 


distribution, we also have 


Sereilary 3.2 


co 


max P = 1- i f {l - p(u,v) }°£, (u,v) du dav 3 2c) 
E°) me 


=O 


subject to 


fo ep (u,v) -audy = A; 
(3.9) 
Uw) Zz in 
where 
2 2 
Ee ne 
1 20% 20% 
f, (u,v) = eae (3.10) 
7,9, 
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has the solution 


2 2 
IL u V 
1 -exp{-o + —— (== + —5)} (u,v) ¢« E 
2 (n-L) on 92 
mo, Vv) = (3.22 } 
0 Gig) cE. 
where 
E = {(u,v)|u*/o% + v*/o2 § 2(n-1)9) 
@ is positive and is determined bv 
2? A/2n(n-1)o, 0, = 0. 12) 


The maximum value of the salvo kill probability is given by 


Pee en es Pll (mal) (t ae 9h) (emia 


Proof 
The equation (3.11) is readily obtained from (3.4) in 


Meoposition 3.1 with £, (u,v) Given wove oe 0) .and C,; 


—i__ g in-lie (3.14) 


Beosctitucing (3.11) in (3.6) and changing variables into 


fe,0) by 
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= y2(n-1) or cos 8) We Y2(n-1l) or sin 9 (e315) 


we have 


ie 
r 


J f p*(u,v) duav 
E 


1 2 
20 (n=!) G..0 il (Das Zr ar 
eV 0 


an(n-1)o, 0, (p-l+e °) 


Piocerenecing (3.10) and (3.14) in (3.7) with (3.15), 


vo 
A -(n-1) > 27 2 
B= =——e- +5————(n-l)o co f -(n-1)r“__-(n-1) 6 
270, 0. 2700, uwvi'y (e e zee Chae 
= l- ee a + (n-1)¢ - A/2mo_ a_ } 
EL 
Midis ftOrmula together with (3.12) gives (3.13). C 


Bor simplicity, let us call the function p*(u,v) the 
(Walsh) optimal p function, and the associated salvo kill 
probability P* the Walsh optimal P. It should be noted that 
the constraint in the revised problem is weaker than that of 
the original problem, and so the Walsh optimal P is an upper 
bound to the maximum salvo kill probability of the original 
problem. It should also be noted that the Walsh optimal p 
function is not feasible in the original problem (being non 


meee On Only a finite set) 1f£ dispersion is present. 





B. APPROXIMATE FORMULA FOR THE OPTIMAL BALLISTIC DISPERSION 


Walsh developed the theory stated above, but he did not 





extend his theory beyond the results (3.4)-(3.7). The impor- 
tant part left by Walsh and subsequent workers is to bridge 
the gap between the revised problem and the original prob- 
lem. Walsh got the optimal p function as given in (3.4), 
whereas we need the optimal density function f(x,y). The 
material in this section bridges that gap. 


Let us define three double Fourier transforms as follows: 


iz oo © i (w,utw.,v) 
p*(wi,W5) = i i aa rake (aL I USubs cht 


a 00 0 i(w,x+woy) 
Dw, ,W5) = f f e Die 7a; (326 ) 


5 00 <0 i (w,X+w>y) 
f, (WW) = f J e £, (x,y) dx dy 


Then from (1.1) we have 
2, P a zs 
p (Wy 1W>) D(w, ,W5) £., (wo) ; CSL) 


In the Walsh theory, D(x,y) is a given function, and 
the Walsh optimal p function is given by Proposition 3.l. 
Therefore, both D(w, ,Wo) and p*(w, ,0,) in (3.17) are now 
known functions. Thus the Fourier transform of the "optimal" 


ex, Y) can be calculated by 
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~ p* (w/w) 
‘co. (Saks) 
2 iv“e BY ) 

GD 


e * Y) e s 
and we will get £, (x,y) aS its inverse Fourier transform. 


This might have been the end of Walsh theory, but there are 
difficulties. First, we do not have any guarantee that 

£5 (x,y) 1s a density function, because Walsh solved the 
restricted optimization problem subject to p(u,v) 2 0 but not 
£. (x,y) 2 0. Further, even if this problem were solved, there 


would still remain the problem of feasibility of the f. func- 
tion. Walsh suggests in his paper that it would be possible 
to obtain a real-world dispersion function £5 (x,y) as the 
theory requires, at least approximately, by careful design 

of the ammunition, but this may be asking too much of the 
Manufacturing process. We assume below that the function 

£5 (x,y) is normal, with only the variance oS subject to 
Memuteacturing control. 

The problem, then, is to determine or approximate the 
optimal ballistic dispersion ooh given the damage function 
Miepy) and the function p*(x,y) from (say) Proposition 3.1. 
In the following, both D(x,y) and p*(x,y) are assumed to be 


even functions of x and y respectively, and to have "moments" 


of all degrees. The lethal area is denoted as A. 


Me xD(x,y)axay = { ff yD(x,y)dxdy = --> =0, (3.19) 


om OO =~ OD = CO = OO 
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ff xp* (x,y) dxdy = ff yp*(x,y)dxdy = +++ = 0, 
f f D(x,y) dxdy = A, (3.20) 
fi fi Eee cca ek = 1,2, 
fj f_ x *Kp*(x,y) dxdy < ©, k = 1,2, 
etc. 


For the sake of generality, we actually consider the non- 
symmetric case where the optimal dispersion may be differ- 


ent in the x and y directions. Let the two dispersions be 
2 Z 





eeeand o . Then 
x oy 
2 2 
ae on 
20% 20 
£, (x,y) ea ; Cae) 
x y 
with the double Fourier transform 
52 2 62 Z 
el See 
@., (iy, ot) 2 (3522) 
pee |. = 2 ° 
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Expanding the exponential function in the double Fourier 


transform of D(x,y), and integrating term by term, we have 


x 00 oo i (w,x+w.y) 
D(w, 14.) = | ji e . Zs Bia ny ) cxdy 


= L =; f i (w, xtwoy) ID (x,y) dx dy 


oO [o.e) 


5 x jar x“D (x,y) dx dy 


( 
> 
‘ai 

} 

| 


ie @) oO 


iL 2 
m _ J y°D(x,y) dxdy + 2.0] 
where the third line 1s obtained because D is even. There- 


fore, for sufficiently small values of ws and ws, 


D(w, Ww) = Dy (W 1 1W>) (3523) 
where 
2 
Dp (wy 15) = Aexp [--> ZF (ee eo, ed dy 
wy 2 co fore) 
= = i { y7D (x,y) dx dy] . (3.24) 
Similarly, 
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olka (W1W>) = Po (W5 Wo) (25) 


fo. small we and ws, where 


2 
~ GU) CO CO 

be oh 2 
Po (Wy 1H) = Aexp [- Ten \ if men (xs, y) dx dy 


= OO = CO 





— f i y-p* (x,y) dx dy] : (2 6) 





Now we determine two parameters of and oF ng IL Moki 


the relation 


~ Bean, Oe 
CRU) Pik ees ey Qo 
2 641 1s ae 

qs 


Then from (3.22), (3.24), and (3.26) we get an approximate 


Benmmula for the optimal ballistic dispersion: 


co co 


2 1 2 + 
= Ix -_ .. ey ery (yy) cada 
(3%62:8)) 
2 i 7 ae. ce 
me iS om i Van or (oo) i= Dear ya) xc | 
where 
( a Lipo g 0 
i) = (3.29) 
| 0 otherwise 
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Note 1 
In this proposition, the approximate formula for the 


Optimal dispersion is derived from the condition that 


£,(w, 10>) coincides with 2A (Ope for infinitesimally small 


wt and ws. It would be better if we could derive the approxi- 


mate formula from the condition that £, Cw >) is nearly 


seas 


* 
equal to £,(w,,4,) Over a wider domain of wt and WS which 


include (0,0). In the next chapter we will give an idea of 


ted 


improving the approximate formula taking this into account. 


Note 2 
Generally, approximation (3.23) with (3.24) is accurate 
only up to the second power of Wy and Wo - If the damage 


function is of Gaussian type, 


Zz 2 
2a° Nac 
Dix V3 = e Y : 
then 
Do (wy 15) = Dw, ,W.) 


is an exact relationship. 
As a corollary, we have the following approximate formula 


for the model with normal errors and Gaussian damage function. 


Consider a salvo firing of n weapons. The density func- 


tions of the random bias and the round-to-round ballistic 


S.1k 





error are, respectively, 


Z 2 
ee 
1 20° 20° 
£, (urv) = 1G co e ' (S202) 
Das 
and 
JL 2 
a ee 
1 20, AIS 
f(x,y) = Ee a ' (ele S Si) 
SENS 
the damage function is given by 
_(u-x)? | (vey)? 
20° 20° 
D(u-x,v-y) = e : (3 234)) 


Then approximate formulae for the optimal ballistic dispersion 


for fixed n and are 
’ Owe iL Or Oy 


(4/0) = [Z(n-1) “9 (0 ,/a,) “(4 /a,) = (n-1) (o,/a,) “1,017 
(37235) 
2A 22 3 ee 7 ee 
er i oe) Ne pe 7%) 710d 


where » is determined by the equation 


aus -l+t+o- ee ae = 0. (3706) 


a2 








An upper bound to the maximum salvo kill probability is given 


by 
ee a or) (lesen) . (een) 


The solution to the Walsh problem with the above men- 
tioned assumptions is already given in Corollary 3.2. By 


use of (3.11), we have 





= f f u*p* (u,v) du dv 
3 , ay a Tear x? 
oj ae ¥2(n-L) oe soy 
- = an f =x-{l-e 2(n-1)) ear 
0 
1059, Ze 
= amet 2) ) = (n-1)A/to oJ 
and 
1 , x*D(x,y) dxdy = a 
]\ a r¥ x 


Substituting these relations into (3.28), we get the approxi- 


mate formula for a (3.35). Similar arguments follow for 


2 


O.. 
ey 


= 





IV. SIMPLE SALVO MODEL REVISITED 


In this chapter we apply the general approximate formula 
for the optimal ballistic dispersion obtained in the last 
chapter to the Simple Salvo Model dealt with in Chapter II 
and investigate its accuracy. 

The assumptions in this chapter are the same as those in 
Chapter II. We consider a salvo of n weapons against a 
target. The impact point error is composed of two parts, 
the random bias and the round-to-round error. The position 
Gf the target (U,V) has a circular normal distribution centered 


at (0,0) with variance ae The impact points of the ith 


ls 
weapon (X,Y) are assumed independent and identically dis- 
tributed circular normal random variables centered at (0,0) 


; 2 
with variance On. Finally, we assume the circular Gaussian 


damage function with parameter : 








7 Ri 
20 
is —a 1 (4.1) 
+ Deon 
il: 
a Xone 
l 205 
tere = S P (452 } 
y ona 
2 
Z 
mG Ser.9, +(v-y) 
Bei 
D(u-x,v-y) = e ; (4.3) 
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items. 5D), Let 


u V ie 
ope = ae = o., (4.4) 
a = a, = a4 


then we have the following: 
An approximate formula for the optimal ballistic disper- 


Sion of the Simple Salvo Model is given by 


a3 /a* = Races =i tae ky 0]* (4.5) 


where 


N 
ll 


(n-1)0%/a* = I1/(e?-149) . (4.6) 


An upper bound to the maximum salvo kill probability is given 
by 
meer 7 fi } (n=) (toe ")} . (4.7) 
Ot ee Ze 
The true optimal O5/a is a function of n and 0,/a eo she 
according to (4.5) and (4.6) the approximate oS /a* is deter- 
Mined by z = (n-1) 04/a*, @ being a parameter. To obtain an 


2 


explicit formula for oS /o. in z instead of the parametric 


form, we have 
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Lemma 4.1 


The series solution to the equation in 6, 
a.) a ob - Wee =. 96 (4.8) 
1S given by 
p? (ae? TR eee ep renee (4.9) 
Proof 
Let the left-hand-side of (4.8) be expanded in 6. 


Then 


noel = 3) ae AO Le - GO ee iaibsn) 


& 
I 


If the terms in the right-hand-side of (4.9) are expanded 


in ¢, we have 


Mey 2 a= 62760 


u = + 
w?/3 = $°/3 - 69/6 + 97/18 - 
722 = S/o) Sama 
i760 = Pea = 

meom which (4.9) is obvious. 0 
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Using the Lemma and the relationship (4.6), we can express 
a3 /a" in (4.5) as a power series in V2/z. If we keep only 
the first three terms in the expansion, we get an alterna- 
tive approximate formula for the optimal ballistic disper- 


sion which we will denote as we 


a : 
¥Z/2 - 1 30 2) ae (4.10) 


es) TS) 
ron re 


a /a = 
with 
Z = (n-1) 0% /a* 


In Table 6 two values of the approximate optimal ballis- 
tic dispersion are compared. One is calculated with (4.5) 
and (4.6), the other with (4.10). Agreement seems quite 
satisfactory, differences appearing at most in the fourth 


decimal digit. The formula (4.10) needs more terms to be a 


Table 6. Comparison of Two Formulae 

7 Ze Ze Daeee 

Z= (n 1)o,/a o4/0 oo/o 
800.0 hes ly) sees 018 OF 
400.0 34 ays) 7 tL eyes) Tal 
2007.0 Sooo Das S0 7 
150050 Sireisiois, 3.8854 
50.0 225069 e067 
Z5 7.0 1. 535 AS ceil 
Ge 0.8475 0.8467 


aS oa" jnsme@anculatea by (4.5) ;,) (4.60), and 


g*/a* is calculated by (4.10). 


a) 





good approximation for smaller values of z, but later we 
feeesee that the case with z less than 25 is not important 


from the viewpoint of application. 


Note 


For n = 2, the formula (4.10) gives 


) a oF oa Vv 2 
om Ae! es o5/a ae 30 a/ o> : 2s lay 


Mimenaotcer it, the exact formula of the optimal ballistic 


dispersion was derived for n = 2. The series expansion of 
2 2 s 
o5/a is 
Zee S ji 


It is interesting to note that every term in (4.11) is smaller 
than the corresponding term in (4.12): The approximate 


value given by (4.11) for n = 2 is always smaller than the 


true value. 


For comparison's sake, we will also consider the classi- 


cal approximate formula by Merritt and King [Ref. 2] who 


approximate the optimal ballistic dispersion by oo where 
2,2 .1-- 
d/o =F Vom (4.13) 
_ SY’ 
Zo = no,/a (4.14) 
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A. ACCURACY INVESTIGATION 

In this section we will investigate the accuracy of both 
our approximate formula and the classical one. 

Table 7A for z = 400, Table 7B for z = 100, and Table 
7C for 2 = 25 illustrate the calculation carried out for the 
accuracy study. For a fixed z, the approximate value for 
the optimal ballistic dispersion divided by ac is a constant 
which is given next to the value of z in the table. On the 
other hand, the true optimal value increases with n: The 
first column is n and in the second column is given the corres- 
ponding true optimal o3/a", the calculation of which is based 
On the necessary condition of optimality (2.13), and requires 
fairly long calculation time. The relative error of our 
approximate formula Ace = (9% - 05) /o5 x 1008 is given in the 
3rd column. Its value is found to always be negative. The 
salvo kill probability associated with o, = 0 is denoted as 


2 


Py eae is an the fourth column. The £f1£th column is the upper 
bound P of the maximum salvo kill probability given by (4.6) 
and (4.7). In the sixth column, we have the true maximum 
salvo kill probability designated simply as P which is calcu- 
Maeea by (2.21), (2.22) with the optimum oS given in the 

2nd column. The entries in the 7th column are the salvo kill 
probabilities associated with the approximately optimal dis- 
persion of The notation Pp. is used for it. The last column 


is the relative loss in the kill probability which would be 


caused if we were to use the approximate value of instead of 


Se) 





QuQunQquwnaquw vg 


NOPAOHDNAAL OO 





Table 7 Accuracy of Approximation (4.196) 
A Z= 400, oaJ/«? = 8.60 

On" /oxt NOE To Ee 2 Pa 

9.34 -7/.9 9.0225 6.0414 6.9419 6.84190 
9.46 ~9.1 8.8634 9.1484 8.1456 8.1454 
9.57 -19.2 0.1988 0.2777 0.27/59 0.27084 
9.76 -11.9 0.2824 9.5229 9.5981 6.5867 
9.84 -~12.6 0.2484 8.6230 0.6956 8.6938 
a 9 2 -13.3 Be2os! 0.7861 0.6873 9.6851 

Z = , Olax? = 3.89 

ou /x? Bice ap Ee P P Ee 

4.41 -11.8 9.1763 8.2947 0.2872 8.2866 
4.48 -13.4 OG. 2828 8.4952 g.4809 6.4794 
4.55 -14.6 0.3877 9.6558 9.6371 6.63508 
4.61 -~15.7 8.4817 Oed727 6.7529 8.7563 
4.66 -~16.6 0.5641 0.8534 9.8348 9.8321 
4.72 -17.4 Gi O352 6.9071 9.8919 6.8885 

2 Gelar = 1253 

Oo /a2 Ai Ora wy Pe P P Sees 

1.77 -13.6 8.0573 8.8668 6.9652 4.9651 
ees 2 Sg = Oo 2a 7 Beas oe 8.2663 0.2657 
36 -17.8 8.3706 8.4822 09.4674 9.4658 
1.94 -19.3 8.5876 0.6497 Deane? 8.6274 
Lopes: -22.4 0.7812 8.9885 8.8912 g.8884 
i oo 23.1 9.8356 8.9435 9.9295 8.9271 
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the true optimal ve AP = ease 5 PLOUs sf it 25 noted that 
in the case with a large z value the salvo kill probability 

is much improved by employing the optimal ballistic disper- 

Sion. For instance, when z = 400, the ratio P/P 4 1s around 

2.5. But the improvement is not so much when z gets small. 

Meee Case with z = 25, the ratio is at most 1.25. 


As to our approximation, we observe the following. 


(1) The approximate of is always smaller than the true 
value, at least for the values of z investigated. 

(2) The discrepancy between the approximate and the 
true values’ is larger for smaller z. 

(3) For a fixed z, the optimal ballistic dispersion of 
increases as n is increased. Thus the discrepancy 
grows when n increases for a given z. 

(4) However, we may say that the discrepancy is not so 
serious because the kill probability associated with 
the approximate optimal ballistic dispersion is not 
so different from the maximum value. The relative 
loss given in the last column is at most 0.4% which 


occurs when z is 25, a less important case. 


mEwiicewo Ls presented here to iltustrate the properties 
l and 2. The ordinate is the relative error Ao®, lonbye alike 
needs several words. The relative error do* varies with n 
for a fixed z as seen in Table 7. To get a representative 
Ao* value for a z, ann is picked for the given z such that 


the maximum salvo kill probability is close to 0.6, an 
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arbitrarily chosen figure. Then the corresponding Ao* 1s 
chosen for that z value, for instance -12.63 for z = 400, 
and similar procedure follows for other z's. The bottom 
curve in Figure 5 shows this relationship. Later we will 
discuss other curves in the same figure. 

From the figure we see that the relative error is always 
negative and its absolute value is decreasing in z. The 
z-~dependence is, however, relatively small, and we learn 
its bias to the negative side is a characteristic of this 
approximate formula. 

Similar investigation was carried out on the classical 
approximate formula (4.13). Table 8A for ti 400, Table 8B 
for Z = 100, and Table 8C for Z = 25 correspond to Table 7, 
the only difference being that ZO 1s tabulated, rather than 
Zz. From the tables we cbserve the following. 

(1) The approximate oe /a* is larger than the true value 
except for the cases of very large z,'S- 

(2) The discrepancy between the approximate and the true 
value is larger for smaller Zo: This tendency is 
Srmilar tO that Of Our approximation. 

(3) For a fixed Zs the optimal ballistic dispersion of 
increases as n is increased. Thus, the discrepancy 
gets smaller when n is increased for a given Z less 
than 800. 

(4) However, the discrepancy is not serious because the 


kill probability associated with the approximate 
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Optimal dispersion is not so different from the true 
maximum. The relative loss is at most 1%, and is 


very small when oS /a is large. 


The top curve in Figure 5 is the relative error do* of 
the classical approximate formula, where the abscissa should 
be read as Zo in this case. In spite of its simple form, it 
is indeed an excellent approximation, in particular for large 
Z,'S- The Z dependence of the relative error Ao® is, how- 


ever, more sharp than ours. The middle curve in Figure 5 will 


be discussed in the next section. 


B. IMPROVEMENT OF THE APPROXIMATION 

In the preceding section, we investigated the accuracy 
of the approximate formula given by (4.10), and it seemed 
that (4.10) gave a lower bound to the true optimal ballistic 
dispersion. The reason it gives a possible lower bound, and 
a method of obtaining improved approximations will be con- 


sidered now. 


* 


, and others in (3.16) are functions 


The Walsh optimal f 


or xX, 
es ne eee) AS) 
1 2 
for the Simple Salvo MOdel, and so let us denote them as 
ee, (=) instead of £5 (wz 1W5), SECU WLtemewEentce notation, (3.18) 


is 
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£5(x) = p(x) /D(x) . (4.16) 


Equation (3.27) by which the approximate a /a* is determined 


reads 
fi{x) = e = Po lx)/D) (x) , ae) 
and 
D(x) = Dy (x) ; (amelie) 


because the damage function is assumed to be Gaussian. 
The function p*(x) is the double Fourier transform of 
the Walsh optimal p function given in Corollary 3.2 with 


oe oy - OL and Pg (x) is an approximate formula for it 


—~ 


given by (3.26). Figure 6 illustrates what the curve of p*(x) 
is like. The curve decreases aS x 1S increased from zero, 
changes its sign, and then swings back to positive with small 
amplitude, and so forth. On the other hand, the approxi- 
mating function Py (x) is equal to p* (x) at x = Q, has the 
same tangent there, and decreases exponentially. Py (x) is 
larger than p* (x) for a wide interval of x as is seen in 
Figure 6. 


~e 


Py (x) > 1. Come 0 < X < Xp 
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Mmeretore, from (4.16) and (4.17), we have 


Pe 
“x0 /a - 
é > £. (x) ; Ox = ex 


which suggests that of foesMa wor euaan) ne Leven optama l 
value. 
One of the simplest ways to get an approximate formula 


would be to define og? 


b 
-8x0"/a* 2 s 
e - Po (x) /Dy (x) (4.19) 
Mmem Oo < 6 < l, or 
2 2 
Oa = oR 7c ; OF aeons (4.20) 


and to determine the parameter 9 empirically to achieve 
satisfactory agreement between the approximate oe and the 


true a Note that the definition of of corresponds to 


2 b 


~ tad 


approximating the p*(x) function by Py (x) with an exponent 
1/8 times the exponent of Po (x) - 
Bor example, let 1/6 be 1.15. Then, 


2 
— ee ee Eo oO, ° a PN) 


In Table 9A, Table 9B, and Table 9C, we present the accuracy 


investigation of (4.21). Accuracy is found much improved. 
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The maximum absolute relative error Aor reduces to 12%, and 
the relative loss in the salvo kill probability is, at most 
only 0.07%. The middle curve in Figure 5 is the relative 


error Ao. 
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V. SALVO MODEL WITH NON GAUSSIAN DAMAGE FUNCTION 


In this chapter we deal with a salvo kill model with 
circular normal errors and circularly symmetric general 
damage function. Under this general assumption, calculation 
of the salvo kill probability is very difficult, and time- 
consuming computation is needed for it. Here we propose an 
approximate method for calculating the salvo kill proba- 
bility. By this method computation is very easy and yet 
fairly high accuracy is expected. 

Assumptions in this chapter are as follows: A salvo of 
mn weapons is fired. The impact point error 1s composed of 
two parts, the random bias common to all n weapons and the 
round-to-round error. We adopt a Cartesian coordinate sys- 
tem (x,y) the origin of which coincides with the center of 
impact of the n weapons. The random position of the target 
is denoted as (U,V) with respect to this coordinate system, 
and (U,V) is assumed to be circular normal centered at (0,0) 


with variance of: 


2 
uo +Vv 
20 
e . (Se) 
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£, (u,v) = 


The impact points of the ith weapon (X,Y) are independent 


and identically distributed circular normal random variables 


gk 





centered at (0,0) with variance of: 
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The conditional kill probability of the target by a single 
weapon given that the target is at (u,v) and the weapon impact 
point is (x,y) is called as the damage function D. [In this 
chapter, D is assumed to be a function of the miss distance 
r= Ete + (voy) - only: 


ee ae (5.3) 


Further, the lethal area is denoted as A, and it is assumed 


that there exist finite "moments" of all degree as defined 


meme ( 5.5) , 
Bey Dix) rdr = A , (5.4) 
0 
E(x?) = ca f rd Di@e) Gene EO , Isao aes. «©6294 5) 
0 


The conditional kill probability by a single weapon 


given (U,V) being (u,v) is given by (1.1), repeated here as 


p(u,v) = f f D(u-x,v-y) f. (x,y) ax dy : (556) 
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peeetehne salvo kill probability (1.2) is renamed as (5.7). 


CO CO 


n 
me= LF -f |; {1-ptu,v)} f, (u,v) dudv . (557) 


—_OO — 


In the Simple Salvo Model in Chapter II with Gaussian 
assumption on the damage function, we have a simple expres- 


mom tOor p(u,v) : 


7 ut+yv 
42 2(a* +04) 
p(u,v) = —+ ze : (25 6) 
a” +O 


2 


In the general model here, p(u,v) given by (5.6) is a com- 
plicated function in Vut sv‘, auGscOmoucatien ©£ the Kill 
Siagoadility is difficult. To overcome this difficulty, 
we propose an approximate formula for p(u,v), and then 
present a recursive method for computing the salvo kill 


probability with this approximate p(u,v) function. 


A. APPROXIMATE FORMULA FOR p(u,v) 

In this section, an approximate formula for the condi- 
tional kill probability p(u,v) is presented. As a corollary, 
weewill get an approximation to the so-called Circular 
Coverage Function. 


Define the double Fourier transform of p(u,v) as 


x co vo i(w,utwov) 
P(w, Ws) == if f e Etuny jeau dy | (5.9) 
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with similar notations being used for D(x,y) and f(x,y). 


Then, from (5.6), we have 
Tne Fourier transform of the circular normal density (5.2) 
is 

£, (w) 14>) = e ‘ (S20) 


In the Fourier transform of the damage function, let the 


exponential function be expanded and integrated term by term, 


a 00 co i(w,xtwoy) 
D(w, 7W>) me if f e D(x,y) dx dy 


= y if if — fe (wx +woy) °D(x,y) dx dy 


0 =00 —<o a) 


eas yrmmetry, the odd terms are zero, so 


(w,xtwoy) ““D(x,y) ax dy 


D(w, 1w.) 1 
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Therefore, 


‘ = (-1)* (2k- mR J o 2%) (yay 2) * 
D(w,1W.5) /A Ll + = ions OM Eaie 15 ( ) (wy +w 5) 
(Se 
since 
(2k-1) (2k-3)--+-1 nS 
5 2k(2k-2)--+-2 : a 
l 7m 2k 
5s J (w, cos@+u, sine)“"de = 
0 
1 k = 
Mier rtirst few terms in (5.11) are 
a ne 2 2 2 li 
D(w,,w5)/A = 1 a Bir ) (wy) +w,) + ee ) (wf + W - 
Now let us define 
‘ pete ~c(wi+w5) 
D, (w,/4>) — 7(w, tw5jie (5 


and determine two constants b and c such that (5.11) and 
foe 2z) COlncide up to the fourth power of ws 'S. Then we 


get 


te ey oe (5 
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fret o>.10) and (5.13) are substituted in (5.9), the result is 


2m 
S Z 2 
~ -——(W_+wW_) 
- - b Dy 2 2 1 2 
where 
s¢ = of - b+ Ee (515) 


meee b given by (5.13). 


Next, the inverse transform of pi fw).w ) is calculated: 
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oo co ~i (uw, +vw,) ~ 
Eau, v) = 5 f f @ P, (W)1W5) dw, dw, : 
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The inverse transform of the first term in (5.14) 
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is easily calculated as 


_ utsy* 
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mine secona term in (5.14) is 
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Zee 2 
p, (u,v) = a {1 = (1-47 te - CSL 7) 
27S S PS. 


= is noted that (5.17) 1s exact when the damage function 
1s Gaussian with parameter a; (5.17) coincides with the exact 
formula (2.6) because in this case A = Amer, b = QO, and 


eee 2 
S = O5+a < 


The cookie cutter damage function differs substantially 
from the Gaussian, so consideration of the former provides 
a good test of the robustness of (5.17). The cookie cutter 


damage function is 
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D(r) = (S23) 
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otherwise 
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The probability 


eu,v) = jf f D(u-x,v-y) £,(x,y) dx dy 
_(x-u) “+ (y-v) * 
ae 
1 2 
= 770. j S : e Gx dy 
x +y Sa 


is called the Circular Coverage Function. In the literature, 
the notation P(a/o,,r/o,) is used for it, where r = Wee awe, 
A standard method of calculating the Circular Coverage Func- 
tion is the algorithm given by Brennan and Reed [Ref. 10] 
based on an infinite series expansion, but if the accuracy 
required is not too high, it is convenient to use closed form 
approximations. Standard handbooks [Refs. 11,12] recommend 


the next formula for a/o. ALY 


_ Hi Dy 
P, (a/o5,X/05) = aoe : (on 9) 


with 
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Poe of ey ae (5.20) 


Peceraing to (5.17), the Circular Coverage Function 


_(xer) *+y? 
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Pla/o,,r/o,.) =s-—- ff e 2 ax ay (5.21) 
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can be approximated by 


2 
es 
2 2 => 
eayo,,t/o.) = Ss{1~ S(1-=,)te 7 , (5,02) 
25 S 25 
where 
mee. 2°/4/3 ; s* = of ta (1 -1/¥3)/4 (5.23) 
Since A = Ta“, ES) = aA PD, and ne) = Boe for the 


cookie cutter damage function. 


mierabile 10, the Circular Coverage Functio P(R,d) is 
meeraced for R= 0.1 (0.1) 2.0 and d = 90.0 (0.5) 2.5. In 
a (R,q) = P, (R,d) - P(R,d) is shown. 


For small R, the error is very small; for instance the error 


Table ll, the error A 


is less than nome when R £ 1.0. But when R is increased 
beyond 1, the error grows rapidly. 

In Table 12, the error A, (R,d) = P_ (Rd) - P(R,d) of 
enblg approximate formula is given. The accuracy is found much 
better than the above-mentioned formula. The absolute error 


4 < 2 


for d = 0 is less than 5x10 ° when R = 1.0, less than 10 


even when R = 1.6 and usually shows much higher accuracy. 


B. RECURRENCE FORMULA FOR THE SALVO KILL PROBABILITY 
If the conditional kill probability p(u,v) by a single 


weapon is given in the form (5.17), we can easily compute the 
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-6.0992 9.8096 
-9.0993 9.8909 
-G.90897 9.0BdB 
-8.09012 9.8080 
-2.8G19 9.9280 
-8.90031 -G.9090 
-2.8848 -9.9901 
-9.8072 -9.8203 
-9.90104 -9.8807 


-0.9199 


-9.8821 


Error A, (R.d) 


se eh) 


-0.G00G 
-9.8000 
-9.880d 
9.8009 
9.9200 


9.8090 
0.GB000 
0.82800 
9.9000 
9.90091 


9.8001 
0.8803 
9.98085 
9.0898 
9.09812 


9.8019 
@.6827 
9.8939 
9.9853 
9.9071 


2.90 


-9. 8800 
-8. 90002 
-G.998G 
-8.90800 
0. BBGO 


9.BBOo 
6.8800 
9.0900 
0.8002 
9.8BBo 


9.9999 
0.9991 
8.9982 
0.8803 
9.99886 


9.0BB9 
0.8815 
9.0823 
0.092834 
9.0949 


Ogaey | 


-2. 9800 
-9.0009 
-9 . 9900 
-G. 8802 


--9. 8000 


-2. 282888 
-8.8802 
-2. 388 
-3.G2BBRB 
-8.080G 


-9.5990 
-8.G9088 
-8.909091 
-2.90901 
-2.9001 


-8.8G02 
-0.99002 
-9.989082 
-9.9001 

9.99081 


salvo kill probability in a recursive way. We first observe 


Piece Ene Salvo kill probability 


co co 


P= l1-f f {1 -p(u,v) }£, (u,v) duav 


= CO = 


where ff, (u,v) is the circular normal density (5.1) and 


Zz 
ee 
2 2 
piu,v) = (urB ay, a 2S : r= BC aye ; (S28) 
AS 
Can be expressed as 
P = J1-= K(n,0) , USNS 745) 
where 
mtn, 0) = J {l - Ma peye 2 hea s dtc, (5.26) 
0 


and where o s*/o?. This is simply a matter of substitution 


and the introduction of polar coordinates. We then have 


Eeoposition 5.1 


tae function 


co 


K(n,j) = f tl - (ut Bt/pje 70 pPen I/PtL) Cae cee 
0 


Be 





which has two integer arguments and three parameters u, 8, 


and 9p satisfies the following recurrence relation: 





K(n,3) = a K(n-1,4) + apap acu)? : EET: (lie, analy 
nee. See = 0 ply. rome) 
with 
K(0,j) = = nice ae (5.29) 
proot 


Baaem (5.27), 


co 


K(n,j) = f {L- (ut Bt/p)e S/ Oy Rte SIZ PHIE ay 
0 


Lo @) 


el om (ue Bt/pye /Pyn hy, st Bt/pye e/a” I/P+1) Cay 
0 


The first term in the right-hand-side is K(n-1,j). The second 


terms is, when integrated by parts, 
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oO 


- f {1- (u+Be/pye 7 PHP Ty + Be/pye OMe” F/P +2) Cay 
0 


CO 


Se 1 - (ut Bt/o)e “9 }Re I/PtLI Et, 
0 


= =(3/p+l) f {l-(u+ Re Peon 2) Ose Ey/oreil) te ate 
0 


oO 


Semen rey eno pte 7 %— C7 e ay 
0 


= 2£¢1-y)" - S22 K(n,3) - BK(n-1,5+1) 
19 
moieretore 
K(n,j) = K(n-1,3) + E(i-y)" - J x(n,j) - @k(n-1,5+1) 


from which (5.28) is immediate. From the definition of 


Mee), 1¢ 2S Clear that (5.29) holds. 


Note l 

We use the notations uw and p here in the same way as 
in Chapter II without any subscripts because no confusion 
is expected. Their definitions are different, but they are 


just the counterparts of yw and 9 defined in Chapter II. 


Note 2 
It is worthwhile to note that the recurrence relation 


feze tas a Similar form to (2.17). The p function (5.24) 
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involves an extra parameter 8 which was missing in the Simple 
Salvo Model in Chapter II, and therefore, it is necessary 


to introduce another index j. 
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VibweCOOKEE, CULTER SALVO MODEL 


Consider that n weapons are fired at a target in a 
salvo. The impact point error is composed of two parts, 
the random bais and the dispersion. We adopt a Cartesian 
coordinate system (x,y) the origin of which coincides with 
the center of impact of the n weapons. The random position 
of the target is denoted as (U,V) with respect to this 
coordinate system, and (U,V) is assumed to be circular normal 


centered at (0,0) with variance of. 





2 
. e 1 ; (Ga 


ib 





eo UL Ae) 
1 Zi 


The impact points of the ith weapon (X.,¥,) are indepen- 
dent and identically distributed circular normal random 


variables centered at (0,0) with variance a5 





Du 
ee yan. 
D 
7 205 
£,(x,y) = 7 (6.2) 
2710 5 


The kill probability by a single weapon conditional on 
mre = (u,v) and (X,Y) = (x,y) Is a function of u-x and 


v-y, D(u-x,v-y), and in this chapter it is assumed that 
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I tee (u-x) * + (v-y)* s — 


Sie} 


Biu-x,v~-y) 


0 otherwise , 


with lethal area 


In the following, we call this model the Cookie Cutter 
Salvo Model. 

The conditional probability that the target is destroyed 
by a single weapon given the random position of the target 


being (u,v) is 
p(u,v) = f { DE THI) a (07) CRG : 


ema with (6.2) and (6.3), we have 


_(u-x) 7+ (v-y) 


2 
Zo 
eou,V) = = 2 - dxdy = P(a/o,,Vurtv*/o,) , 
Pde Se PA 


UO oer +y <a 





(6.4) 


where P(R,d) is the Circular Coverage Function. 


The salvo kill probability 


oO 


2 fe f _{1-p(u,v) }°£) (u,v) du dv , 
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with (6.1) and (6.4) is readily rewritten as 


oo 


S1= 
= 1 - fee eee, 120; /0,))"e shee (6.5) 


We carried out a number of computations of (6.5) with 
various sets of parameters Oi: Oar a, and n by the Simpson 
feevea. A robust algorithm for computing the Circular 
Coverage Function P(R,d) for wide ranges of R and d was needed. 
For this purpose, we used the Brennan-Reed formula in a modi- 


mea form. 


Brennan-Reed Formula: 


f 


2 
a 7(2n) ’ ine = IL 


nw 
~ 
— 
ae 
2 
! 
ae 
lI 


2 
Shay. 
kp ee ——- 
= Uae yop a R/2 jn n= 1,2 
Sn Gye aad 
2 
Gg = 1 - eos ve 
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Modified Formula: 


De, es 2 

Mee) = e° /* Jy (1-e 74% jh , (6.7) 
nN ig! 
iO 
Jaa = 6 fam ASL. 
af “peg 124 vee 

i ln 
Co = Co 1 + ane ae Sine es saat (62-3) 
a /a = R“/(2n) = IL, 2 
n’ n-l ‘ 6 
Co = ao ae RL a 


The modified formula is especially of use in preventing 
underflow possible in cases with large R and d. 

In their review paper, Eckler and Burr presented a short 
table of the salvo kill probability of the Cookie Cutter 
Salvo Model. But the accuracy is, as they admitted, sup- 
posedly low because it is obtained by Monte Carlo simula- 
tion. In Table 13, we give the correct values. The figures 
in parentheses are from the table by Eckler and Burr. 

In the following sections, comparison of the Cookie Cutter 
Salvo Model and the Simple Salvo Model is dealt with and 


then the optimal ballistic dispersion is investigated, where 
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Table 13 


ald, 


0.4 


0-8 


A/a, 


0.4 


e228. 
(8.2089) 
g.458 
(0.447) 
9.605 
(g.595) 
Q.715 
(0.786) 
9.801 
(8.792) 
9.865 
(9.856) 
G.912 
(0.988) 


Calo, 


0.50 


9.228 
(9.235) 
8.587 
(9.582) 
0.836 
(9.832) 
9.952 
(8.951) 
G.990 
(8.988) 


0. /0, 


0.75 
Ole2s 


GERI) 


0.551 
kOe 55 le) 
O7772 
(0-773) 
0.894 
(9.892) 
Oa 
(9.956) 
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O77 
(8.187) 
09.514 
(doe) 
0.764 
(O76 3 ) 
025796 
(8.893) 


(0.957) 


(8.958) 


(Or. 97 5) 


Salvo Kill Probability : Kisi/(Eckler-Burr) 








the approximate formula presented in Chapter III and the 
approximate method of computing the salvo kill probability 


developed in Chapter V are found to be useful. 


A. COOKIE CUTTER VS. SIMPLE SALVO MODEL 
The Cookie Cutter Salvo Model differs from the Simple 
Salvo Model of Chapter II only in the damage function. [In 


Chapter II, we assumed 


_ (arx) 24 (v-y) 7 


; 2 
D(u-xX,v-y) = e = - 


with lethal area pe For comparison's sake, the lethal 


area of the two salvo models are now set equal 
a nt; an 6.29) 


The two models share the common parameters Jy, 9 Cow ancl mm, 


>! 
only being different with respect to the shape of the damage 
function. The quantities associated with the Cookie Cutter 
Salvo Model will be given suffix c in the following when 
necessary. 

It is interesting to compare the salvo kill probabilities 


of the two models. We begin with three extreme cases in which 


the salvo kill probability has closed form expressions. 


(1) The case n=l. 
The formula (6.5) with (6.4) is easily integrated and 
give 


o) 02 





P = Jl-e : ore) 


The corresponding single shot kill probability under the 


assumption of the Gaussian damage function is 


P = ee C672) 


Proposition 6.1 
The Cookie Cutter Salvo Model with n = 1 gives a higher 


Eeeate Shot kill probability than that of the Simple Salvo 


Model. 
P ee ee CEA IZ) 
Cc 
Proof 
=< 2 2 2 : 
Recall that e = Jel sc) fOr xeo vO .- Let. 4 /(o, +95) aes 
mato .l0) and (6.11), then (6.12) is obvious. U 


Generally, in the salvo with n > 1, this inequality 


does not necessarily hold as will be shown in the following. 


(2) The case o, = Q 


Let nT > 0 in (6.5), then we have 


Zee 
-noa [05 
P = Jl-e : (6.03) 
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The corresponding formula in the Simple Salvo Model is given 


by (2.23), namely 


2 = We Wao) 1G A) 
3) The case J5 = Q 
Let 0, > O in (6.4). We get 
-a*/o8 
P = —- re ; Useless) 


The corresponding formula in the Simple Salvo Model is given 


mei 2.24): 
af _ got) al 
B= 2s (SN GSE Saar woe28) 
with 
Le a= a*/o% 


Mmeeposition 6.2 
Between the salvo kill probabilities of the Cookie 
Cutter Salvo Model and the Simple Salvo Model with the same 


parameters, the following relationships hold: 


ee one = O, then 


P Se ae (6.17) 
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> Ze. «2 
im < -P uyseie J /0 - L/ho ; ‘fe cdkey) 
where Xo is a positive solution to the equation 
n 
h- j} An(l + A/j) = 0. (6.19) 
j=l 


Proof 
The proof of (6.17) is just the same as that of Proposi- 


tion 6.1. As to the second assertion, let 


ee eI!) 


mg@eeco.16). Then 


lel nN 
joy = j») IAn(l + A/j) , do/ame— a) 1/7 (+A) 
Wek en 
2 2 
peice g{(0) = 0, dg/dal > 1, a°g/dr~ < 0, we have 
=0 
Te) enn rete \ Sak 


where ho is given by (6.19). The relation (6.18) is 


immediate. 0 


Shs. 





As to the general case, we must compare the two kill 
probabilities numerically. To illustrate the general idea, 
Table 14 is given. In Table 14a, Po for n = 4 is tabulated, 
and the corresponding Simple Salvo Model kill probability 
P is in Table 14B. It is observed that P's in the column 
a4 /o* =u ware Larger than P's as stated in Proposition 6.2. 
P's in the first row corresponding to oS /a* = 0 is smaller 
than P's except for the case of /a* = 0, since 1/5 = 0.216 
bor n= 4. 

As a general tendency, P is larger than P when o, is 


2 


relatively larger than oa Es 1s smaller than P when o 


2 


as is seen in the upper 


1° 
1s considerably smaller than Ors 
right corner of the table. 
PeeeeOrPLIMAL BALLISTIC DISPERSION 
We carried out a golden search calculation to get the 
optimal ratio a3 /a* for a number of cases with different 
combinations of n and o*/a°. The results are tabulated in 
Table 15. Table 15A gives the optimal 05/0 of the Cookie 
Cutter Salvo Model, and Table 15B gives the associated 
salvo kill probability, namely the maximum salvo kill proba- 
bility for a given set of n and of /ae. It is noted that this 
calculation is very time-consuming. The corresponding 
values for the Simple Salvo Model are tabulated in Table 16. 


Before starting comparison of the optimal ballistic 


Gispersion of the Cookie Cutter Salvo Model with that of the 
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Ox /¢2 


PmWWNNdYRFeE QQ 


QoaqnQquwna neg 


Tx 1K? 


Q 
Q 


mWWNYHEe FQ 
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Table 14 


9.2 


1.9882 
1.899 
9.982 
9.931 
8.865 


9.736 
0.681 
9.632 


0.0 
809 


-938 
.879 
892 
749 
-684 
-634 
-599 


QA QAAAAAre 


Salvo Kill Probability, n = 4 
A Cookie Cutter Salvo Model, Pe. _ 
OY M2 

2.0 4.9 6. 8. 
2) NS 8.221 9.154 8.118 
8.685 oj aleys) 8.285 8.224 
8.626 9.425 8.320 9.256 
8.615 @.435 8.335 Bie272 
Geos? 8.432 4) 4 Sei) 9-279 
9.565 8.423 0.337, O23 
Ge538 8.412 G23933 Oe27 2 
G5 1) Bee eis) 8.326 8.276 
@.486 8.385 BSR cule, Baz 72 

B Simple Salvo Model, P. 

oa? 

2.2 4.0 6. 8.2 
8.594 9.382 9.289 Baz 1 
8.621 9.421 Orreriey, O.2253 
8.613 9.431 0.331 8.268 
8.589 G.427 8.334 8.274 
8.569 @.417 B33) 8.274 
8.531 9.404 9.325 8.271 
9.593 9.399 mo. ol / G.267 
8.476 W375 9.399 9.262 
@.452 8.361 9.398 9.256 
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18.0 


- 895 
-185 
-214 
229 
555) 7 
-249 
- 241 
- 240 
2a 





Table 15 Optimal Ballistic Dispersion and 


Maximum Probability ( Cookie Cutter Salvo Model ) 


A Optimal Ratio Ore ay oe 
y\ 
qf x? 2 4 g 16 aD 64 
2 g.4 eo) We 7. 27 al ee dk 
4 g.8 1.6 2.6 3.9 5.8 8.6 
6 5 2.1 Be 4.9 yee a 18.5 
8 0 Daa) 3.8 5.6 8.3 ae) 
10 1.8 2.9 44 6.4 9.3 eae 4 
12 Zak a3 4.9 7.2 O.1 We 7 
B Maximum Salvo Kill Probability 
Pe Max 
on ee 2 4 8 16 32 64 
2 G.476 G.626 Q.782 G.9985 g.973 a.996 
4 G.300 G.435 G.601 G.769 g.899 ¢.971 
6 G.224 8.339 Q9.492 Q.667 9.826 r.935 
8 g.182 9.281 9.420 9.599 9.762 °.897 
19 g.154 G9.241 9.368 9.531 9.708 © .859 
12 9.135 9.213 9.329 G.484 9.661 f.823 
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Table 16 


Optimal Ballistic Dispersion and 


Maxmum Probability ( Simple Salvo Model ) 


A 


oJ? 2 
2 0.5 
4 Oia 
6 G5 
8 Oey 
19 9.9 
12 1.3 
B 
ee rot 2 
2 9.450 
4 9.293 
6 9.217 
8 9.173 
19 6.144 
12 9.124 


Optimal Ratio 62 /«?. 


NNFF Fe @ 
WF O — F OD 


6.622 
9.431 
9.334 
0.274 
BO. 233 
§.204 


PmWWND 
QQ aArF OVO W 


P max 


9.780 
Bae o ede, 
6.489 
9.416 
9.363 
69.323 


ee) 


16 


NM UI & W NY 
W~] QW Bb 


16 


9.964 
9.768 
6.665 
9.588 
9.528 
9.481 


iO © ~} OV UI Ww 
« ee e @ @ GJ) 
UNA WO NY 


Maximum Salvo Kill Probability 


a2 


5) SS 
§.899 
9.826 
9.761 
6.7087 


ON RAeys) =) 


64 


5.6 
eae 
18.9 
ae ee) 
Wee, 
14.1 


8.996 
Oe o7 
D223 
2.896 
8.858 
Gee Phe) 





Simple Salvo Model, we recall the approximate formula (3.22) 


given in Chapter III. 


co oo 


z ee J x? {p* (x,y) = 1B) (ea) cls dy ,0]” Goro) 


In the case of the Cookie Cutter Salvo Model, 


i <a IB 5 We) eb tai ao pZ ; 


Pie 
— 


whereas in the Simple Salvo Model, 


f x? Dix, c-ecy. -— a, 


=O =O 


Pir 
—,, 


Therefore, (6.20) states that the optimal ballistic dispersion 


of the Cookie Cutter Salvo Model is larger than that of the 
Simple Salvo Model approximately by 0/2: 
2 2 25 ae 
05/2 - o5/a = cee (Gree) 
Now we subtract entities in Table 16A from those in Table 


15A and get the difference of /a* - 05 /a* which is tabulated in 


Table 17A. The values are positive with only one exception, 


and the predicted value 0.5 is observed in the upper right 
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Table 


17 


Difference of Optimal Ratio 


0.926 
0.097 
9.007 
0.009 
0.9190 
9.011 


and Simple Salvo Model 


9.004 
0.004 
9.095 
0.097 
9.898 
0.9089 


Pemax 


0.982 
9.002 
0.003 
0.904 
0.085 
0.006 


OR 


' Difference of Maximum P'‘s 


P max 


Comparison of Cookie Cutter Salvo Model 


Cac /a* - 

16 32 
B.5 B.5 
B.5 BES 
B.6 B.5 
Z.6 B.6 
B.7 B.6 
B.7 B.6 

16 20 
G.9G1 Z.QBO 
G.GH1 B .2OO 
G.GH2 D. 200 
G .0G2 G.0G1 
G.GB3 Z.0G1 
G .0B3 B. B02 


oR ac ee 


2.080 
2 .8OD 
°.901 
£.991 
g.9O1 
C.991 





Semmer Of the table. In the below left corner of the table, 
the Figures are not close to the value 0.5, but these cases 
Gerrespond to low salvo kill probability. 

In Table 17B, the difference of the maximum probabilities 
of each model is tabulated. The figures are all positive, 
and lead to the speculation that P > P in all cases 


Cmax max 


despite the fact that Ee may be smaller than P when o.4 is 
not chosen optimally. The author has not yet succeeded in 
obtaining a proof of this speculation, but it might be the 


Beemmeerpart for n > 1 of Proposition 6.1 for n= 1. 


In Chapter V, we developed an approximate method of 
computing the salvo kill probability which can be applied 
to the Cookie Cutter Salvo Model. The salvo kill proba- 
mercy 1S given by 


P = 1 - K(n,0) , (S125) 


where 
ogy = ff fl - (i+et/p)e 7p] (I70TIE a (5.20) 
0 


satisfies the recurrence relation, 





x n 2 : 0 ee ng _ ; 
K(n,j) = seepes iS “ares Hy) nrjtpe ‘2 aie co) 
mo = 1.2) ° 3 a ee 





with 


tO) j=) ey = ao). (5.29) 


ieempaLameters are given by (5.22), (5.23), and (5.24). 


(a“/s*) (1 - a*/2v3 s*) , 


u = 
g = 0 JAG 54 , (6.22) 
s¢ = a - oo Gi - 1/73) /2 


UE 


Using these formulae together with the golden section 
search calculation, we obtained Table 18. Table 18A gives 
the approximate optimal ballistic dispersion of the Cookie 
Cutter Salvo Model, and the associated salvo kill probability 
jmeein Table 18B. 

In Table 19A, the difference between the optimal of Jae 
given in Table 15A and the approximate om in Table 18A 
is tabulated. The difference decreases towards the upper 
mean. corner of the table. 

In Table 19B, the salvo kill probability of the Cookie 
Cutter Salvo Model associated with the approximately optimal 
ballistic dispersion mentioned above is given. It corres- 
Bends to the maximum P given in Table 15B. In spite of the 
discrepancies observed in Table 19A, the approximate of /a 


gives almost the same kill probability as the true ono gives. 


OS 





Table 18 Optimal Ballistic Dispersion 
and the Associated Salvo Kill Probability 


( Cookie Cutter Salvo Model, Approximation ) 


A Splimal RatiowGn. 76... 

Ti la® 2 4 g 16 32 64 

2 9.5 1.0 eee | Pe | 4.1 Sal 

4 Q.7 ea) Zeon 3.9 5-6 8.6 

6 1. 1.9 Sie 4.8 had EOS 

8 ie erz ae 3.6 S65 8.2 2750 

19 hee! 2.2 4.1 6.2 Ore 13.4 

12 4 ee 2.0 4.5 6.8 G.O 14.7 

B Maximum Salvo Kill Probability 
Pa 

oi / ok? 2 4 8 16 32 4 
2 G@.475 9.625 9.782 8.995 9.973 O22 l6 
4 9.297 9.432 9.609 9./68 9.899 6.971 
6 9.218 9.334 9.499 8.665 9.826 0.935 
8 8.173 9.274 9.416 9.588 G9./62 8.396 
19 9.144 8.234 8.363 8.528 G9.707 9.358 
ji 9.124 9.204 9.323 9.481 9.659 fe) sl 7a 
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Table 


Error 


19 Accuracy of the Approximation 


( Cookie Cutter Salvo Model ) 


Poeeeie (OPC Tincaencaw toe 


QaAaq gVgand 
MW BWNrF @ 


QAaqaQqaqagQd 
H&WNOr FH 


2 
Coe [«x* ars oe la? e 


—~ 


QAaAaRgAaB 
NONrFrFAAQ A 


a2 


QAAQQQ 
Meer Qa @ 


B Salvo Kill Probability Associated with 


@.475 
0.382 
Be223 
O.181 
Oia 15 3 
9.134 


0.626 
8.435 
9.339 
0.280 
G.241 
9.212 


0.782 
0.601 


9.492 | 


9.428 
8.368 
8.328 
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16 


0.905 
0.769 
0.667 
9.596 
9.531 
0.484 


the Approximately Optimal Ballistic Dispersion. 


PC ) 


32 


eos 
0.899 
0.826 
0.762 
8.788 
9.661 


€4 


G.996 
BASF sh 
Oa935 
Gea27 
Geeta he, 
O-o25 





Viena vO MOOrIS WITHOUT CLRCULAR SYMMETRY 


The approximate formula for the optimal ballistic dis- 
persion presented in Chapter III covers the salvo models 
without circular symmetry. In this chapter the formulae 
Pomecomputing the salvo kill probability of these models 
are dealt with. 

Suppose that a salvo of n weapons is fired against a 
maeget. The impact point error is composed of two parts, 
one being common to all the weapons, and another being 
round-to-round dispersion. The random target position (U,V) 
with respect to the coordinate system (x,y) is assumed to 
be elliptical normal centered at (0,0) with variances 
of and a” respectively. 

2 
1 20° 20 
£, (u,v) = +e : (751) 


2 
Toe, 


/< 


dh 


The impact points of the ith weapon (X.,Y.) are indepen- 


dent and identically distributed elliptical normal random 


variables centered at (0,0) with variances of and Ove 


respectively. 


2 
= 


qs 
20, 


, 


Ze 


KE NO 


eS) 0) gi area = (7.2) 
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The kill probability D by a single weapon conditional on the 
target being at (u,v) and the weapon impact point being (x,y) 
is a function of u-x and v-y only, and is given by the 


elliptical Gaussian function with parameters Ot. and By 


Specifically, 
D D 
Uo) Gay) 
2 2 
20 Zo 
D(u-x,v-y) = e a Y (G75.3°) 


T€ might be an appropriate model, e.g., for a ground target 
VS. a weapon impacting the ground obliquely, because in this 
case the scattering of the fragments and thus the conditional 
kill probability D is by no means circular symmetric as 
assumed in the previous chapters. For brevity, we call 
this model the Elliptical Normal Salvo Model. 

The kill probability of a single weapon conditional on 


(U,V) being (u,v) is given by 


p(u,v) = f f D(u-x,v-y) £, (x,y) dxdy ; 


menusing (7.2) and (7.3) we get 


2 2 
- u Z V 
2 ne 
2(a° +07) 2(a +o) 
a. 
payy) = ————" _____ — Y ¥ 4) 
|) ee Zane 
Oar oi) MGR) 
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The salvo kill probability is then given by 


2 2 
_Uu _ Vv 
1 ~ = n 208 20% 

ee >) | t1-plu,v)} e dudv . = (7.5) 


oo 00 
From this equation, Grubbs derived the following formula: 


Proposition 7.1 (Grubbs) 
Tieektil probability of the Elliptical Normal Salvo 


Model stated above is given by 


N 
Jel ny, 73 PxPy 
-1 * ar Paar Rea TE _ SE Ty fe} 
2, ae: \ 35,7 G5) (eG) 


Pp = 
where 
és oye: 2 2 2 
u = a0 / hog or Gh) Mak ’ 
C75) 
ey ee SS Oe ae 
ae Fo ey Sg Ce 


The Grubbs formula is not suited for calculation when n 
is large, since it is an alternating series. To overcome this 
difficulty, Breaux and Mohler [Ref. 13] gave a method for calcu- 
lating the kill probability based on an expansion of (1-z)” 
in Jacobi polynomials rather than as a binomial series. The 
series is found to converge to the true value with less than 
n terms, which is very attractive for the calculation of 


salvo kill probability with a large n. 
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The author has tried to find alternative ways of calcu- 
lating the salvo kill probability which are effective even 
for the cases with large n's, but at present, the obtained 
results awe not promising. The following propositions are 
given not as efficient algorithms for calculating the salvo 
kill probability, but as possible hints for developing more 
Smerci1ent ways of calculation. 

For the Elliptical Normal Salvo Model, we have a 
proposition similar to Proposition 2.5. We omit the proof, 


since it closely parallels that of Proposition 2.5. 


BeOpOSition 7.2 
The salvo miss probability Q(n) of the Elliptical Normal 
Semvyo Model characterized by (7.1), (7.2), and (7.3) with 


the weapon number n, 





yoo. 8 pOesesp.t 
Qin) = —22¥ ff (1-ve S * )"e * Y asat (7.8) 
satisfies the relationship 
ay n N= kek 
oye) (CL) Wee (7.9) 
k=0 
where 
2 2 
/ Ok DCeReCe apes =P.,¢ 
q(k) = SEY yf (a-eS "Pe * ¥ asat . (7-10) 
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Pimetewcase Otesimple Salvo Model, q(k) had a simple 
recurrence relationship, but we have not succeeded in 


finding such a relationship in this case. 


Corresponding to the recurrence formula for the salvo 
kill probability presented in Chapter II, we have the next 
propositions: The series given in Proposition 7.3 is an 


infinite series, but for a small §/p it will converge rapidly. 


EeeposSition 7.3 
The salvo miss probability of the Elliptical Normal 


Salvo Model with n weapons is given by 


sts ! 
Qin) = V1i-8*/p* y SA'S (570)** gin, 2k) , (7.11) 
oO) ee 
where 
S(n,k) = a f Gamer! Oy = e © * ai 7 
a- ~ 20 
with 
be ae ae oe 
ie xX x 
cgi tee 
Cu oy 
(eZ) 
2 2 2 2 
(eas ee Q +090 
§ = = = 5 re +) ' 
oa oy 


JIG) 





and 


Cielo 
x Y 


iD 
(oy + oF) (a, + By) 
Proof 
Mimnv ec) let Ss = ¢cosds, t = rsinsd, and (pe) 72 = 0, 
(9,,-0,)/2 = 6, Then 


OSes On G = 2 ly ds 8 BSE 29) 


and eherefore O(n) is 


yee: = ane Bee ee ee 
on) = —— f (niet ace tae {f eameco S ae}x Gs 
0 0 


ren) 2 2 
2/0 Py : Glcate = yoenos 1) (6x°) mor a; 


where I, (2) is the modified Bessel function of the first 


kind with order 0, and has an expansion formula 


fe) 


ie 
Rae) = L Car) 


k=0 


(5) ae (Te): 


Thus 


1G AL 





tI 


Q(n) /o,0,/0 f d-ue*/°)"e*z, (st/p) dt 
Y 0 





22:98 i 


2k a =-t/o,n -t, 2k 
age (6/p) “" sae f (l-ue pyre vt 
ONY = 2% (k1) 7 (ic yal 0 


=v ait 


U 
Proposition 7.4 


The integral J(n,k) which has two integer arguments and 


two positive parameters u and o, 


Bimyk) = gf (l-ne 7°)" @ > &* at (7.14) 
0 


satisfies the following recurrence relation. 


ee = 2s = 
finm,kK) = es JI(n-1,k) + as J (apk=l)s Cis) 
= 2 eee eta? Me eel revs 
a ee = ee eo =. 
J(n,0) a n+o J(n Hi.) a n+p (1 W) ’ n I 2 vein, eet (756) 


UO hs), i = Cieley) 


Pr2 





Proof 


mime) = =f (i-ywe “°)"~ e@ ~ t* at 
= a f eee ems heor tomer ae 


0 


The first term in the right-hand-side is J(n-l1,k). If 


the second term is integrated by parts, then 


(© @ 


J(n,k) = J(n-1,k) - = 2 Glemen! yews t*| 
: 0 
ey el a he a 
eee 0 
6) -t/fo,n.-t _,k 
me oo a be )-e Gc 7ee 


Therefore, for k # 0, we have 

S(n,k) = J(n-1,k) + £ a(n,k-1) - £ o(n,k) 

From this relation (7.15) is immediate. If k = 0, 
3(n,0) = gJ(n-1,0) +& (1-y)" - 2a(n,0) , 

and therefore, (7.16) results. The relation (7.17) is obvious. 


U 


iis 





It is noted that all the terms in the expansion (7.11) 
are positive, and that the integral J(n,0) is the salvo miss 
probability of the Simple Salvo Model dealt with in Chapter 


II. Therefore we have 


Corollary 7.5 
A lower bound to the salvo miss probability of the 
Elliptic Normal Salvo Model characterized by (7.1), (7.2) 


and (7.3) with weapon number n is given by 


Qn(n) = V1 = 8*/p* Q(n) , (7.18) 


where Q(n) is the miss probability of the Simple Salvo Model 


by n weapons with 


2 2 Z 

ioe eS ov es 

8] 2) 

u V 
+O to° 
5 eee wae sy, 

on? 52 52 ’ 
u V 
oe 


a g) Zz Z 
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